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i r*| Abstract: In this paper we continue to pursue a goal of finding an effective theory for high energy interaction in 
QCD based on the colour dipole approach, for which the BFKL Pomeron Calculus gives a low energy limit. The 
£N| key problem, that we try to solve in this paper is the probabilistic interpretation of the BFKL Pomeron Calculus 
^ in terms of the colourless dipoles and their interactions. We demonstrate that the BFKL Pomeron Calculus has 
two equivalent descriptions : (i) one is the generating functional which gives a clear probabilistic interpretation of 
y—( the processes of high energy scattering and also provides a Hamiltonian-likc description of the system of interacting 
C~ J dipoles; (ii) the second is the Langevin equation with a specific noise term which is rather complicated. We found 
■^j" that at high energies this Langevin equation can be reduced to the Langevin equation for directed percolation in 
the momentum space if the impact parameter is large, namely, b » I/fc, where k is the transverse momentum 
of a dipole. Unfortunately, this simplified form of Langevin equation is not applicable for summation of Pomeron 
^ loops, where one integrates over all possible values of impact parameter. We show that the BFKL Pomeron calculus 
. with two vertices ( splitting P — > P + P and merging P + P — > P of Pomerons ) can be interpreted as a system 
bf colourless dipoles with two processes: the decay of one dipole into two and the merging of two dipoles into one 
dipole. However, a number of assumptions we have to make on the way to simplify the noise term in the Langevin 
equation and/or to apply the probabilistic interpretation, therefore, we can consider both of these approaches in the 
present form only as the QCD motivated models. 
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1. Introduction 

The simplest approach that we can propose for high energy interaction is based pi ^ on the BFKL Pomeron || 
and reggeon-like diagram technique for the BFKL Pomeron interactions |[ [| 0]. This technique, which is a 
generalization of Gribov Reggeon Calculus ||, can be written in the elegant form of the functional integral (see || 
and the next section). It is a challenge to solve this theory in QCD finding the high energy asymptotic behaviour. 
However, even this simple approach has not been solved during three decades of attempts by the high energy 
community. This failure stimulates a search for deeper understanding of physics which is behind the BFKL Pomeron 
Calculus. On the other hand, it has been known for three decades that Gribov Reggeon Calculus has intrinsic 
difficulties |J that are related to the overlapping of Pomerons. Indeed, due to this overlapping we have no hope 
that the Gribov Reggeon Calculus could be correct in describing the ultra high energy asymptotic behaviour of the 
amplitude. The way out of these difficulties we see in searching for a new approach which will coincide with the 
BFKL Pomeron Calculus at high, but not very high, energies (our correspondence principle) but it will be different 
in the region of ultra high energies. In a spirit of the parton approach we believe that this effective theory should be 
based on the interaction of 'wee' partons. We consider, as an important step in this direction, the observation that 



has been made at the end of the Reggeon era 110, 11 [12|, that the Reggeon Calculus can be reduced to the Markov 



-i 



process [jl3[ for the probability of finding a given number of Pomerons at fixed rapidity Y . Such an interpretation, 
if it would be reasonable in QCD, can be useful, since it allows us to use powerful methods of statistical physics in 
our search of the solution. 

The logic and scheme of our approach looks as follows. The first step is the Leading Log (1/x) Approximation 
(LLA) of perturbative QCD in which we sum all contributions of the order of (as hi(l/x)) n . In the LLA we consider 
such high energies that 

a s hx(l/x) w 1; while a s < 1 (1.1) 

It is well known that the LLA approach generates the BFKL Pomeron (see [[| and the next section) which leads to 
the power-like increase of the scattering amplitude (A oc ^ n =o, v =o) with u(n = 0, v = 0) cx as). 

The second step is the BFKL Pomeron Calculus in which we sum all contributions of the order of 

a% 1 \ n a% 1 , 
-tto — 1 — n FvT therefore — -, — - -r w 1 (1.2) 

where as = asN c /ir. 

The structure of this approach as well as its parameter has been understood before QCD E3] and was confirmed 
in QCD (see [|l[ ||, ||, |[ ||, [?], |l5[ |l6|). This calculus extends the region of energies from \n(l/x) ss 1/ces of LLA to 
ln(l/x) w (1/as) ln(Nc/as)- The BFKL Pomeron Calculus describes correctly the scattering process in the region 
of energy: 

1 fN 2 \ 1 1 

— In ) « In- « -y 1.3 

as V a s J x a s 

For higher energies the corrections of the order of (a 2 s ln(l/x))" should be taken into account making all calculations 
very complicated. 

Our credo is that we will be able to describe the high energy processes outside of the region of Eq. ([b^) , if we 
could find an effective theory which describes the BFKL Pomeron calculus in the kinematic region given by Eq. ( |l.3| ) , 
but based on the microscopic degrees of freedom and not on the BFKL Pomeron. In so doing, we hope that we can 
avoid all intrinsic difficulties of the BFKL Pomeron calculus and build an approximation that will be in an agreement 
with all general theorems like the Froissart bounds and so on. Solving this theory, we can create a basis for moving 
forward considering all corrections to this theory due to higher orders in as contributions, running QCD coupling 
and others. 

The goal of this paper is to consider the key problem: the probabilistic interpretation of the BFKL Pomeron 
Calculus based on the idea that colourless dipoles are the correct degrees of freedom in high energy QCD j^]. We 
believe that colourless dipoles and their interaction will lead to a future theory at high energies which will have the 



BFKL Pomeron Calculus as the low energy limit (see Eq. (1.3)) and which will allow us to avoid all difficulties of 
dealing with BFKL Pomerons at ultra high energies. 

Colourless dipoles play two different roles in our approach. First, they are partons ('wee' partons) for the BFKL 
Pomeron. This role is not related to the large N c approximation and, in principle, we can always speak about 
probability to find a definite number of dipoles instead of defining the probability to find a number of the BFKL 
Pomerons. The second role of the colour dipoles is that at high energies we can interpret the vertices of Pomeron 
merging and splitting in terms of probability for two dipoles to merge into one dipole and of probability for decay of 
one dipole into two ones. It was shown in [|~7j that P — > 2P splitting can be described as the process of the dipole 
decay into two dipoles. However, the relation between the Pomeron merging (2P — > P) and the process of merging 
of two dipoles into one dipole is not so obvious and it will be discussed here. 

This paper is a next step in our programme of searching the simplest but correct approach to high energy 
scattering in QCD in which we continue the line of thinking presented in |ll| [lj], ^l[ . The outline of the paper 
looks as follows. 



In the next section we will discuss the BFKL Pomeron Calculus in the elegant form of the functional integral, 
suggested by M. Braun about five years ago Q. In the framework of this approach we find a set of equations 
for the amplitude of n-dipole interaction with the target. We show that the recent intensive work on this subject 
[ p2| , |25| , |20| ] confirms the BFKL Pomeron Calculus in spite of the fact that these attempts were based on slightly 
different assumptions. 

In section three we demonstrate statistical interpretation of the theory with interacting Pomerons. The one-to- 
one correspondence between BFKL Pomeron calculus and Langevin theory is found by showing how full Lagrangian 
generates stochastic Langevin equation with a peculiar noise term. In toy model with zero transverse dimensions 
this noise is reduced to one typical for directed percolation. Unfortunately the complexity of the noise term restricts 
the practical use of Langevin equation in this form, and one should look for further simplifications. One of them 
is to assume that impact parameter is much larger than any dipole size in the system (see Eq. ( |3.20 )). Using this 



assumption and going to momentum space we rewrite our theory in the form of Langevin equation with a noise term 
proportional to the field (directed percolation universality class). It should be mentioned that large impact parameter 
approximation is unapplicable for summation of Pomeron loops, where we integrate over all possible values of the 
impact parameter. 

Next, we discuss an approach based on generating functional. We show the equivalence between generating 
functional approach and the BFKL Pomeron calculus in the kinematical region Eq. ( |l.3| ) that leads to a clear 
interpretation of the BFKL Pomeron calculus as an alternative description of the system of interacting colourless 
dipoles. The interrelation between vertices of the Pomeron interactions and the microscopic dipole processes is 
considered. It is instructive to notice that the generating functional approach leads to a feedback to the BFKL 
Pomeron Calculus restricting the integration over Pomeron fields in the functional integral by the range < $ < 

In the fourth section we suggest a practical way of building the Monte Carlo code to solve the equation for the 
generating functional which can be a basis for consideration of the multiparticle production processes. 

In conclusion we are going to compare our approach with other approaches on the market. 
2. The BFKL Pomeron Calculus 

2.1 The general structure of the BFKL Pomeron calculus 

We start with a general structure of the BFKL Pomeron calculus in QCD. The BFKL Pomeron exchange can be 
written in the form (see Fig. [l]-l) 

-2 

A(Fig.@-l) = V u (g)Gp(r 1 ,r 2 ;b\Y 1 -Y 2 )(^V d <x -jL exp (w(n = 0, v = 0) (Y 1 - Y 2 )) (2.1) 

^ c 

with u>(n = 0, v = 0) cx asXi — Y 2 = ln(l/x) and ® denotes the all needed integrations. 

It is easy to understand the main parameters of the BFKL Pomeron calculus by comparing the contributions of 
the first 'fan' diagrams of Fig. [l]-2 with the one BFKL Pomeron exchange. 

This diagram has the following contribution 

A(Fig. 0-2) = / dY'V u (g)Gp(r ll r';b\Y 1 -Y')(g)r(1^2)G 2 P (r' 1 r 2 ;b\Y'-Y 2 )(g)Vi (2.2) 



x Vu , V * y 2 J. exp (2 w(n = 0, v = 0) (Y 1 -Y 2 )) cx (&) exp (2w(n = 0, v = 0) Y x - Y 2 ) 
oj{n = 0, v = 0) 



where r± and r 2 are the sizes of the projectile and target dipoles while {r'} denotes all dipole variables in Pomeron 
splitting and/or merging. 

-2 

One can see that the ratio of this two diagrams is proportional to j& exp (u>(n = 0, v = 0)(Yi — Y 2 )) which is 
the parameter given by Eq. fll.2Q . When this ration is about 1 we need to calculate all diagrams with the Pomeron 
exchange and their interactions (see Fig. [lj-a, - Fig. 0-f ). All vertices, that are shown in Fig. 0, has been calculated 
in ^, ||] and they have the following order in «s 1 : 

-2 - 2 

uj(n = 0,u = o) cx a s ; r(l -» 2) oc ^ ; r(2 - 1) cx ^ ; 

- -2 - - 

T(2^2) cx ^ ; r(2^3) cx ^ ;V d oc ^; (2.3) 

Using Eq. ( |2.3D we can easily estimate the contributions of all diagrams in Fig. 0. Namely, 

A(Fig.|-3) cx L 2 (Y) ; A(Fig. 0-4) cx ^YL 2 (Y) ; (2.4) 

A(Fig.&5) cx j^L 3 (Y); A(Fig. 0-6) cx |u 2 (y); (2.5) 

with L(Y) = («|/iV c 2 ) exp(w(n = 0,v = 0)Y). 

As we have mentioned in the introduction the BFKL calculus sums all diagrams at such a high energy that 
parameter L(Y) is of the order of 1 (see Eq. (|l.2j )). In this kinematic region we need to take into account the 
diagrams of Fig. 0-1 , Fig. 0-2 and Fig. [TJ-3 (see Eq. (pj|),Eq. ( |2~2| ) and Eq. Indeed, diagrams of Fig. 0-4 and 

Fig. 0-6 ( see Eq. (jD|) and Eq. tfHfy) are small since (a|/A^ c 2 ) Y < l/N 2 < 1 in the kinematic region of Eq. (|Di|), 
while the diagrams of Fig. 0-5 (see Eq. (pH[)) are small at L (Y) m 1. 

The first conclusion that we can derive from this analysis that in the kinematic region where L (Y) ~ 1 we 
need to take into account all diagrams with T(l — > 2) and T(2 — > 1) vertices while the diagrams with T(2 — > 2) and 
r(2 — > 3) vertices give small, negligible contributions. 



However, if L (y) cx A^ c one can see from Eq. (2.1) - Eq. ( |2.5| ) that all diagrams give so essential contribu- 
tions that we have to take them into account. Indeed, for such L(Y) A(Fig. |lJ-3) oc N 2 , A(Fig. H-4) oc asY > 
l,A(Fig. 0-5) oc Nl and A(Fig. 0-6) oc a s Y > 1. 

It is interesting to notice that the vertex T(2 — * 2) oc ag/N 2 can be neglected even at such large values of L (Y). 

Finally, we can conclude that the first step of our approach can be summing of the diagrams with T(l — * 2) and 
T(2 — * 1) vertices in the kinematic region L(Y) 1=3 1 or Y 1/as In (iV^/as). 

However, we would like to stress that we need to make an additional assumption inherent for the BFKL Pomeron 
calculus: the multi-gluon states in i-channel of the scattering amplitude lead to smaller contribution at high energies 
than the exchange of the correspondent number of the BFKL Pomerons (see more in pi]]). This statement is 
supported by the fact that numerous attempts to find the intercept of these states being larger than the intercept 
for multi-Pomeron exchanges[E5l have failed. 



2.2 The path integral formulation of the BFKL calculus 

The main ingredient of the BFKL Pomeron calculus is the Green function of the BFKL Pomeron describing the 
propagation of a pair of gluons from rapidity Y 1 and points x[ and x' 2 to rapidity Y and points x\ and X2 2 ■ Since 

1 In Eq. ( ^.3| ) we use the normalizations of these vertices which are originated from calculation of the Feynman diagrams. In the dipole 
approach we use a different normalization (see below section 3 and 4) but all conclusions do not depend on the normalization. 

2 Coordinates Xi here are two dimensional vectors and, strictly speaking, should be denoted by Si or Xj. However, we will use notation 
Xi hoping that it will not cause difficulties in understanding. 



V u V d m->2) r(2->i) r(2->2) r(2->3) 




Figure 1: The BFKL Pomeron interactions and the examples of the diagrams of the BFKL Pomeron calculus in QCD. The 
solid line describes the Pomeron exchange while the double line stands for the dipole. 



the Pomeron does not carry colour in i-channel we can treat initial and final coordinates as coordinates of quark and 
antiquark in a colourless dipole. This Green function is well known p6|, and has a form 

G(x' 1 ,x' 2 ;Y\x 1 ,x 2 ;Y') = Q{Y-Y')x (2.6) 

OO „ 

x dv d 2 x e"( n > v X Y - Y ') \{n,v)E n {x[,x 2 ;xo\v)E*( Xl ,X2;x \v) 

71= — OO 

where vertices E n are given by 

E n (x u x 2 ;xoW) = (^-) h (^K-)' 1 (2-7) 

\X 10 X20j \X W X m J 

where Xik = Xi — Xk, Xi = x^ x + iXi iV 3 ,x* — Xi. x + ixi >y ; h — (1 — n)/2 + iv and h = 1 — h* . The energy levels 
u>(n, v) are the BFKL eigen values 

w(n,i/) = a s (^{\)-Re^(^^+i^\ (2.8) 
where ip{z) = d\n.T(z)/dz and T(z) is the Euler gamma function. Finally 

") = [( n + 1)2 + 4^2] [ (n _ 1)2 + 4jy 2] ( 2 ' 9 ) 

The interaction between Pomerons is depicted in Fig. ^| and described by the triple Pomeron vertex which can 
be written in the coordinate representation || for the following process: two gluons with coordinates x[ and x' 2 at 



3 Xj^ x and Xi tV are components of the two dimensional vector xi on x-axis and y- axis 



Figure 2: The graphic form of the triple Pomeron vertex in the coordinate representation. 



rapidity Y[ decay into two gluon pairs with coordinates x'{ and x 3 at rapidity Y 2 and acg and x' 3 ' at rapidity Y$ due 
to the Pomeron splitting at rapidity Y. It looks as 



7ra| f d 2 X\d 2 x 2 d 2 x 3 



111 

x 12 x 23 x 13 



(G(x;,4;y/|a ;l ,x 2 ;y)^ ll2 )G( Sl ,X3;r|a:i^X3^r 2 0G(x3,x 2 ;y|x3^4^r3 , ) (2.10) 



where 



^1,2 = r \2V\v\ withp 2 = -V z 



(2.11) 



and arrow shows the direction of action of the operator L. For the inverse process of merging of two Pomerons into 
one we have 



2 7ra s f d 2 x\d 2 x 2 d 2 x 3 



N c 



2 2 2 

X* rf^ 1 s-p*-- 
12 x 23 x 13 



G(x'(, 4'; Yl\xi,x 3 ;Y) G(x' 3 \ x'{- Y^\x 3 ,x 2 ;Y) (T l<2 G(x 1 ,x 2 ;Y\x' 1 ,x' 2 ; Y{)) (2.12) 



The theory with the interaction given by Eq. (2.10) and Eq. (2.12) can be written through the functional integral 
as was proposed and developed by Braun in ||. We include a discussion of the basics of this approach for the sake 
of completeness of our presentation. 



Z[$,$+] = J e s with S = S + Sj + S E 



(2.13) 



where Sq describes free Pomerons, Sj corresponds to their mutual interaction while Se relates to the interaction 
with external sources (target and projectile). From Eq. ( 2.10| ) and Eq. ( 2.12| ) it is clear that 



S = J dYdY' d 2 x 1 d 2 x 2 d 2 x[d 2 x 2 $ + (x 1 ,X2]Y)G- 1 (x u X2]Y\x[,x 2 ;Y')$(x' l7 x 2 ;Y / ) 
dy f d 2 x,d 2 x 2 d 2 x 3 {(ri ^ (xi;X2 . y)) $+( XljX3 .Y)*+(x a ,X2;Y) + h.c.} 



Si = 

N c j j u. 12 u^ 2 3 .i.13 

For Se we have local interaction both in rapidity and in coordinates, namely, 



S E = - J dYd 2 xid 2 x 2 {$(xi,x 2 ;Y)T pr (xi,x 2 ;Y) + $ + {xi, x 2 ; Y) T tar (xi, x 2 ; Y)} 



(2.14) 
(2.15) 

(2.16) 



where r pr (r tar ) stands for the projectile and target, respectively. The form of functions r depend on the non- 
perturbative input in our problem and for the case of nucleus target they are given in ||. 



For the case of projectile being a dipole that scatters off a nucleus the scattering amplitude has the form 

Txi,yi;r =T<U(xi,yi;Y) = — — ± — = , 2.17 



where extra as comes from our normalization and we neglect term with r pr in Eq. ( |2.16 ) 



Generally, for the amplitude of interaction of n dipoles at rapidity Y we can write the following expression 4 

The extra factor (—1)™ is due to the fact that in Se the source for both projectile and target, has extra minus 
sign. 

It is useful to introduce the Green function of the BFKL Pomeron that includes the Pomeron loops. This function 
has the form 

G(x 1 ,y 1 ;Y\x 2 ,y 2 ;Y) = / d¥d^ eW>* + l\s E=0 ( } 

For further presentation we need some properties of the BFKL Green function p6[ ]: 

1. Generally, 

G- 1 (x 1 ,x 2 ;Y\x\,x' 2 :Y') = V \ V \ ( A + h) = ( A + H+) V \ V \\ (2.20) 

as f d 2 x 3 x 2 
2ir J x 23 x^ 

2. The initial Green function (Go) is equal to 

2 2 2 2 

G {x 1 ,x 2 ;Y\x[,x' 2 ;Y) = it 2 ln ^ 1 '^ 2 ' In X \ V ^ (2.22) 

■^1,2' ,2 x l,2 ^l'^' 



dY I \dY 



Hf{xi,xr,Y) = / 2 3 2 12 (/(zi.aajY) - f(xi,x 3 ;Y) ~ f(x 3> x 2 ;Y)); (2.21) 
z7r J a; 23 2^3 



This form of Go has been discussed in pfl. In appendix A we demonstrate that this expression for Gq(x\,x 2 ; Y\x[,x' 2 ; Y) 



stems from lo = cu(n = 0,v) term in sum of Eq. (2.6). Only this term is essential at high energies since all other 
terms lead to contributions decreasing with energy. 

3. It should be stressed that 

VlV 2 2 G (x 1 ,x 2 ;Y\x' 1 ,x' 2 ;Y) = (2.23) 

(2tt) 4 {5^{xx-x' 1 )8^{x 2 -x l 2 ) + 6<»(xi- x' 2 )8^\x 2 - x[)) 



4. In the sum of Eq. fl2.6| ) only the term with n — is essential for high energy asymptotic behaviour since all 
u){n, v) with n > 1 are negative and, therefore, lead to contributions that decrease with energy. Taking into account 
only the first term one can see that G is the eigen function of operator L\ 3 , namely 

L 12 G{x 1 ,x 2 ;Y\x' ll x' 2 -Y l ) = TTrT^ G(xi,x 2 ;Y\x[,x 2 ; Y') w G(x u x 2 ; Y\x[, x' 2 ; Y'); (2.24) 



4 Starting from this equation we use notations Xi for the coordinates of quark while yi denote the coordinates of antiquarks. For 
rapidity we will use Y . 



The last equation holds only approximately in the region where v -C 1, but this is the most interesting region 
which is responsible for high energy asymptotic behaviour of the scattering amplitude. 

All properties of the BFKL Pomeron Green function as well as of the functional integral approach to the BFKL 
Pomeron calculus have been discussed (for more information see g, pfifl). 

In the next section we will derive the chain of equations for multi-dipole amplitude in the BFKL Pomeron calculus 
and will show that these equations are the same as ones that have been discussed in framework of dipole approach 
§§§■ 

2.3 The chain of equations for the multi-dipole amplitudes 



Using Eq. ( 2.13 ) and Eq. ( 2.17 ) we can easily obtain the chain equation for multi-dipole amplitude noticing that 



every dipole interacts only with one Pomeron (see Eq. (2.17)) 



These equations follow from the fact that a change of variables does not alter the value of functional integral of 
Eq. (|2.13|) . In particular, = Z[$,$'+] (sec Eq. (|]l|)) where $'+ = $+ + e(x,y) with a small function 

e(x,y). Therefore, 

" D$D<P+e s ^^ = [ D<t>D<P'+e s ^^ (2.25) 



(2.26) 



Substituting $' + = $ + + e(x, y) and expanding this equation to first order in e, we find 



= J D<S>D<S>+e s ^^ x 
J dY dY' d 2 x\ d 2 X2 d 2 x' 1 d 2 x 2 e(xi,X2\ Y) G~ 1 (xi,x 2 ; Y\x' x , x' 2 ; Y') ^{x\, x 2 . Y') 

{(e(x 1 ,x 2 ;Y)I: h2 ) ■ <t>( Xl ,x 3 ;Y)<f>(x 3 ,x 2 ;Y)-{ 



9 9 9 

ry,*. 

x 12 x 23 x 13 



2(p 1 j${xx,xr,Y)) ■ e(x 1 ,x 3 ;Y)^> + (x 37 x 2 ;Y)} - J dY d 2 Xl d 2 x 2 e{x u x 2 ; Y) T tar (x u x 2 ; Y)] 

We redefine the integration variables in the third term as follows 

2(r 1)2 $(n,o; 2 ;y)) • e(x u x 3 ; Y) $ + (a; 3 , x 2 ; Y) —> 2 (Zr li3 $(a:i, x 3 ; Y) ) ■ e(x 1 ,x 2 ;Y)<S> + (x 2 ,x 3 ;Y) 



Using the expression for the Hamiltonian Eq. 2.20| and the Casimir operator Eq. 2.11 we define a new variation 
parameter e(xi, x 2 ; Y)p\ p\. In terms of this variation parameter Eq. 2.26 reads as 



+ p S[<£> ,<£>+] 



(2.27) 



[J dYd 2 Xl d 2 x 2 e{x 1 ,x 2 ;Y)p 2 lP 2 2 ^ A + Hj <S>(x u x 2 ; Y) 

{(e(x 1 ,x 2 ;Y)p 2 p 2 ) ■ $(xx,x 3 ; Y) $(x 3 , x 2 ;Y)+ 



2ira 2 s f nr f d 2 x\d 2 x 2 d 2 x 3 



dY 



o-«2 ,y>2 2 
x 12 x 23 x 13 



2 2 p 2 2 

2(r i3 $(xi,x 3 ;r)) ■ e(x u x 2 ;Y)^^ + (x 3l x 2 ;Y)} - / dY d 2 x 1 d 2 x 2 e(x 1 ,x 2 ;Y)^ T tar (x 1 ,x 2 ;Y)} 

P\P 2 J PiP 2 

We denote the new variation parameter by e(xi,x 2 ;Y) = p 2 p 2 e(xi, x 2 ; Y) and use the property of the initial 



Green function Eq. 2.23 to rewrite -rs in terms of Go as follows 
1 1 P1P2 



-2—2^i,x 2 ;Y) 
PIP 2 



2 {2tt) 4 



G {x u x 2 - Y\x[, x' 2 ; Y)i(x[, x' 2 ; Y)d 2 x[d 2 x' 2 



(2.28) 



Thus, Eq. 2.27 can be written as 



= / D<S>D<S>+e s ^^ 



(2.29) 



j dYd 2 Xl d 2 x 2 e{x u x 2 ;Y) + Hj §{ Xl ,x 2 ;Y) 



2na 2 s 



2(r i Mx 1 ,x 3 ;YY 



dY 



d 2 xi d 2 x 2 d 2 x 3 



"12 ^23 ^13 



{e(xi, x 2 ;Y) <f>(xi,x 3 ;Y) $(x 3 , x 2 ; Y)+ 



2(2tt) 



G (x 1 ,x 2 ;Y\x' 1 ,x 2 ;Y)i{x' 1 ,x' 2 ;Y)d 2 x' l d 2 x' 2 } $+(x 3 , x 2 ; Y)} 



- / dYd 2 Xl d 2 



x 2 \ 2 (2ir) 4 J G o( x i^ x 2;Y\xi,x 2 ;Y)e(x 1 ,x 2 ;Y)d x i d x 2f T tar(xi,x 2 ;Y)] 



Noting that the r.h.s. of Eq. (2.29) should vanish for any possible variation of e{x\, x 2 ; Y) we obtain 



(2.30) 



+ 



2 7ra| 2 



N c (2tt) 4 



2nal 
H - 

N c 

d 2 x[ d 2 x' 2 d 2 x 3 

2 2 2 
•*'l'2" i '2'3 ^VS 



x\ 2 d 2 X 3 
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<f>(xx,x 3 ;Y)$(x 3 ,x 2 ;Y)+ 



(x\, 3 P 2 , P mx' 1 ,x 3 ; Y) ) G (x[,x' 2 ;Y\xi,x 2 ;Y) $ + (x 3 ,x' 2 ; Y)} 



1 



2(2tt) 



d x[d x' 2 G ( X ' 1 ,x 2 ;Y\ Xl ,x 2 ;Y)Ttar(x' 1 ,x' 2 ;Y)} 



where we interchanged (x[, x 2 ) <-> (xi, x 2 ). We notice that the third and last terms are independent of rapidity and 
can be absorbed in the initial condition. This is obvious for the last term which represents the target source. To 
show this for the third term we use the property of the Casimir operator at high energies (n = , v = 0) 

E h3 $(xi,x 3 ;Y) ~ $(xi,x 3 ;Y) 

and the definition of the Green function (see Eq. |2.19| ). This equation will be discussed in the next section in more 
details. We see that the third term results into the product of two initial Green functions which are independent of 
rapidity. 

Now we can use the definition of the amplitude defined in Eq. 2.17 and Eq. 2.1§| to rewrite Eq. 2.30 in a simple 
form 



dTW( Xl ,x 2 ;Y) = as 
dY 2tt 



d 2 z K(x\, x 2l z) 



(2.31) 



(T^\x u z;Y) + T^{z,x 2 ;Y) - T^( Xl ,x 2 ;Y) - (Y; xi,z; z, x 2 ; Y)j 
where kernel K{x, y\z) is defined as 

K(x,y\z) 



(x - yf 



(x — z) 2 (z — y) 2 



(2.32) 



and the Hamiltonian is given by Eq. 2.21 



This equation has a very simple meaning that is clear from Fig. |^. 
Starting from equation 



(2.33) 




Figure 3: The graphic form of equations for the multi dipole amplitude. 

we obtain the equation for the amplitude T^ 2 \ namely, 

dTW{ Xl ,xxx 3 ,x A -,Y) = a^J d 2 zK{xuX2lz) 

(t {2 \xi 1 z;x 3 ,x 4: ;Y) + T (2) (z, x 2 ; x 3 , x t ; Y) - T^ 2) (x 1 , x 2 ; x 3l x 4 ; Y) - 
-T ( - 3 \xi,z;z,X2;x 3 ,x i ;Y)j + 

+ 11? / d 2 x'd 2 X ''r 2 ^ 1 (x u x 2 ;x 3 ,x i \x',x'')Vl,Vl,,T^(x',x''- 1 Yy, 
where function T 2 —,i is equal to 



F 2 ^i(x 1 ,y 1 ;x 2 ,y 2 \x,y) = / d 2 zK(x, y; z) G ( Xl , y ± \x, z) G (x 1 , y ± \z, y) 



(2.35) 



In derivation of Eq. (2.31) and Eq. (2.34) we use Eq. (HH) and Eq. ( Uf) as well as normalization condition (see 
Eq. (2.17)) for the scattering amplitude. These two equations are the same as in This shows that approach 

developed in |22, p7| gives a result similar to that obtained from BFKL Pomeron Calculus. 



Assuming = T<U we obtain the Balitsky-Kovchegov equation |^8|, |2j| . We can do this only if we can 
argue why the Pomeron splitting is more important than the Pomeron merging. For example, this assumption is 
reasonable for scattering of the dipole with the nucleus target. Generally speaking, the splitting and merging have 



the same order in as ( see Eq. ( 2.1C ) and Eq. ( 2.12 ) ). In Eq. ( 2.31 ) and Eq. ( 2.34 ) these two processes look like 
having a different order of magnitude in as, but this fact does not interrelate with any physics and reflects only 
our normalization. However, we will see that for a probabilistic interpretation the correct normalization is very 
important. 
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3. Generating functional and probabilistic interpretation 



3.1 Statistical physics analogy: Langevin equation and directed percolation 



The functional of Eq. ( 2.1 3| ) is reminiscent of the partition function of statistical mechanics. Indeed, the partition 
function has a general form 

Z[H] = e -W F I ff J = f dse -Tk f d X (H(s)-Hs( X )) (31) 



where F[H] is the Hclmholtz free energy. As an example, Eq. ([O]) is written for the system of spins with the energy 
density Tl(s) in the external magnetic field H. The integration is performed over all possible spin configuration in 
the system. 

Comparing Eq. ( pM| ) and Eq. ( 2.13 ) one can see that Eq. ( 2.13| ) describes a statistical system with kT = 1 and 
with F[—T tr ] = —S. The form of Se suggests that H = — T tr plays a role of the external field. 

We demonstrate the idea of the statistical interpretation of our theory by showing how one can generate full 
interaction Lagrangian starting with Langevin equation with a specific choice of a noise term. 

We define a partition function as an integral over all possible field configurations being restricted to satisfy 
Langevin equation with a Gaussian noise 



(3.2) 



S[ I d 2 Xl d 2 x 2 dY{ d ^ Xl d * 2]Y) +H$( Xl ,x 2 ;Y) 



2iras f x 2 2 d 2 x 3 

~NT 



§(xi,xs;Y)${xs,X2;Y) + C(aui , x 2 :Y)} 
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where P(C) stands for Gaussian probability distribution of the noise with < £{xi, x 2 ; Y) >= 0. The noise ((xi, x 2 \ Y) 
has correlation function 



< C(x 1 ,x 2 ;Y)C(x' 1 ,x' 2 ; Y') >= B 5{Y - Y') 5{x x - x[) 6(x 2 - x' 2 ) 
where B is some function of xi, x 2 and Y, and will be determined later. 



(3.3) 



We want to rewrite Eq. (3.2) as a functional integral with a term responsible for merging of Pomerons. To do 
this we introduce a response field $(a;i, x 2 ; Y) and use exponential representation of <5-function 



Z[*,$,C] ~ / P(0 DC D$ exp( ? J d 2 Xl d 2 x 2 dY ^>{ Xl , x 2 ;Y){ d ^ Xl ^ 2 ' ,Y ^ + ffSfa, x 2 ; Y) 
+ 2na l I x^xs Y )<S>( XSiX2] Y) + Cfo, s a ; Y)» 

Performing functional integration over noise C(xi,x 2 ;Y) first, we have 

~ / exp(i / d 2 xid 2 x 2 dY ^(xi,x 2 ;Y){ d ^ Xl '^ 2 ' ,Y ^ + H<S>( Xl ,x 2 ; Y) 



(3.4) 



dY 

277,1 S ' "^p-^xuxs; Y)<S>(x 3 ,x 2 ;Y) + *-B $>(x x ,x 2 ; Y)}) 
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(3.5) 



Here we used properties of Gaussian integral. Now we compare functional integral of Eq. (3.4) with one given by 



Eq. (2.13). It is clear that in order to have proper free and interaction terms we should identify arbitrary <fr(xi, x 2 ; Y) 
and B as 



i 3>(xi,x 2 ;Y)=j% v\ <f> + (xx,x 2 ;Y) 



(3.6) 



-i -i 



and 

The last missing thing is to consider functional integration measure. After the redefinition of the response field the 
measure of functional integration is D{—i v\v\ & + (xi, x 2 ;Y)) instead of desirable D($ + (xi, x 2 ; Y)). But according 
to the definition of the functional integration we have to sum over all possible field configurations. After performing 
change D{— i p\p\ & + (xi, x 2 ; Yj) — > x 2 ; Y)) we still account for all possible configurations up to, maybe, 

an infinite constant, which cancels in any calculations of physical observables. 

Thus, using a Gaussian noise of the form Eq. ( |3.7[ ) we fully reproduce functional Eq. ( 2.13| ), except the terms 
of Lagrangian corresponding to external sources. As it was already mentioned these terms can be absorbed in the 
initial condition of Langevin equation. It should be stressed that correspondence between functional and stochastic 
formulation of the problem was shown on general basis without any additional assumptions. 

In the toy model of zero transverse dimensions the action is given by 



s = J dY f $+(y){ JL - A}$(r) + G 3P ($+(r)$ 2 (y) + $(r)$+ 2 (y)) - $ 



-<P+(Y)T tr 5(Y-Y ) ) (3.8) 

where A and G 3 p are Pomeron intercept and triple Pomeron vertex. 
For this form of action Eq. ( |3.7|) becomes 

B = 2G 3 p$ (3.9) 



leading to correlation function Eq. ( |3.3| ) to be proportional to the field. This means that our theory belongs to direct 
percolation universality class. 

As a final step of our discussion we would like to write Langevin equation used in Eq. ( |3.2[ ) in a slightly different 
way. Because of Gaussian form of probability distribution P(C) we may make the following change of variables 



t(x 1 ,x 2 ;Y)^VBC(x 1 ,x 2 ;Y) (3.10) 
In this case corresponding Langevin equation reads 

d®(xi,x 2 ;Y) 2vra s f x\ 2 d 2 x 3 



dY 



H^( Xll x 2 ;Y)-^- f ^p-^(x 1 ,x 3 ;YMx 3 ,x 2 ;Y)+VB(:(xi,x 2 ;Y)=0 



(3.11) 



and correlation function Eq. (3.3) is given by 

< C(xi,x 2 ;r)C(x' 1 ,a; , 2 ;y') >= 6(Y - Y') 6{x x - x[) S(x 2 - x' 2 ) (3.12) 



In the toy model the Langevin equation for redefined noise Eq. (3.10) is given by 

^p- = A$(y) - G 3P <P 2 (Y) - V2G 3P $(y)CCn (3.13) 

Therefore, the description of the BFKL Pomeron calculus with the help of the Langevin equation for directed 
percolation is equivalent to the functional integral with the action of Eq. ( 2.13 ). In QCD we expect (see |2^, |2l|] ) a 
different form of correlation function ( < £C >oc $(1 — $) ) which belongs to universality class different from that 
of Eq. 0. 



Going back to coordinate dependence we must admit that the form of noise Eq. (3.7) is too much complicated for 
calculations, and we will use Li t2 $>(l, 2) ~ $(1, 2) approximation for our further considerations. We believe that this 
approximation describes well Pomeron dynamics at high energies. The justification of this assumption is given in the 
next section, where we consider Langevin equation for a mixed momentum and impact parameter representation. 



3.2 Langevin equation in mixed representation 

At high energies both n and v are small and one can assume that 



Li, 2 $(l,2)~$(l,2) 



VfV^(l,2)~ — $(1,2) 

X Y1 



(3.14) 



This assumption means that for all BFKL Pomerons in the BFKL Pomeron calculus we can use the diffusion 
approximation, which is equivalent of using Fisher-Kolmogorov-Petrov-Piskounov equation in statistical approach. 
In other words, we assume that all partial energies in the BFKL Pomeron calculus are high. In this case the form of 
the action of Eq. ( 2.13| ) simplifies considerably 



So + Sj = 



d 2 X\d 2 X2 



dY 



J x 12 x 23 x l 3 ^7T 



■ 1 ,x 2 ' ) Y){$(x u x 2 ;Y) - $(x u x 3 ;Y) - ${x 3 ,x 2 ;Y)} 



-^^{^(x x , x 2 ; y)$+(a;i, x 3 ; Y)$+(x 3 , x 2 ;Y) + /i.e.} 



(3.15) 



To justify the assumption of Eq. ( 3.14 ) we want to show how the simplified action of Eq. ( |3 .15 ) generates the 
same BFKL calculus as the full one. To see this we consider the enhanced diagram of Fig. |l|-3 using the full action 
of Eq. (2.13), namely, 



A(Fig.$-3) = J 



d x -^d x^d x^ d x -yd x^d ' 

"J12ZJ12ZJ12 



Y' 



dY' / dY"G{x\,xl;Y\x' l ,x' 2 ;Y')T: v 2 , 



■_, 12 u. 13 ^23 X 1 2 X 1 3 X 23 JO 

xG(x[, x 3 ; Y'\x±, x 3 ; Y")G(x' 3 , x' 2 ; Y'\x 3 , x 2 ; Y")T, V , 2 ,,G{x'[, x' 2 \ Y"\x\, x 2 ; 0) 



Eq. (p.6[) can be written in the following form 



G(x\, x 2 ; Y\x-i, x 2 . Y ) — 



OO r. „ 

£ / dv J d 2 x e^ n ^ Y - Y '^G n {x^xl;x^x' 2 \ V ) 



(3.16) 



(3.17) 



n— — 00 



Therefore 



00 „ „ 

G{x\y 2 -Y\x^x' 2 -Y')*C v . 2 , = J2 J dv j d'xoe^V-^Gnixl^x^x'^Zy^ 



n— — 00 



J dv J d 2 xoe 



u)(n,v)(Y-Y') _ 



n— — 00 



A(n, v) 



G n ( x i: x 2 't x 1t x 2 \ v ) 



dv I d xoe 



w(0,v)(Y-Y') 



A(0,0) 



G n= o(x l , x 2 \ x l5 x 2 \v) 



G(x^, x 2 ] Y\x x , x 2 ; Y ) 



(3.18) 



where we used the fact that Y — Y' , Y' — Y" and Y" — are so large that we can account only for n — term and 
replace A(0, u) by 1 (see Eq. ( |2.9| )). In Appendix [A] we show that replacing the sum over n in Eq. (3.17) by the n = 
term we obtain the Green function of the BFKL Pomeron that satisfies initial condition determined by the Born 
approximation. Using Eq. ( 3.1 8| ) we can rewrite Eq. (3.16) as 



A(Fig.$-3) 



d 2 x' 1 d 2 x' 2 d 2 x' 3 d 2 x'[d 2 x'id 2 x'{ <- Y 



J, 12 X 13 X 23 



r ll 2rr ,ll 2rr tl 2 



Y' 



dY' dY'dx^x^Y^x'^Y') 



X \ 2 X \ 3 X 23 Jo 

xG(x[, 4; Y'\x'{, 4'; Y")G(x' 3 ,x' 2l Y'\xl 4;Y")G(x'{, a#; Y"\x\, 0) (3.19) 
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which is the expression for the diagram of Fig. [j|-3 for the action in Eq. ( 3.15 ) 



Eq. (3.15) can be reduced to a more elegant form in mixed momentum and impact parameter representation. 
We assume that impact parameter is much larger than a dipole size 



Xi + x 2 



3> X\ - x 2 



(3.20) 



Indeed, we expect that the typical size of the dipoles will be of the order of inverse saturation momentum 1/Q s (x), 
while the typical impact parameter of the scattering dipole should be much larger (at least of the order of the size of 
the hadron). Thus, we take b — (xi + x 2 )/2 ~ (x\ + x^)/2 ~ (x 2 + x^)j2. 

We introduce the Fourier transform 



*(xi,x 2 ;y) = x{ 2 \ d 2 k e lk ' X12 $(fc, 6; Y) 



(3.21) 



In terms of fields $(fc, b; Y) and <E> + (fc, b; Y) the action looks as follows 



S = Sq + Si 



\ oY 



2tt 



d 2 k'K(k,k')$(k',b; Y) 



2naj 



{$+(fc, 6; Y)®(k, 6; Y) + $(fc, 6; F)$(fc, 6; Y)} (3.22 



where K(k, k') is the BFKL kernel in the momentum representation, namely, 



d 2 k' K[k,k') &{k',b;Y) 



{k-k') 2 (k' 2 + {k-k') 2 ) 



(3.23) 



{k-k') 2 

In Eq. (3.22) we obtained interaction terms multiplied by a constant vertex. In the derivation of Eq. ( 0.22 ) we used 
condition Eq. ( 3.20 ), in this limit one can neglect momentum transferred along Pomerons resulting into a constant 
vertex. The action Eq. ( [3.2S ) looks similar to the action of the toy model Eq. (3.8). Following the procedure presented 
in Section H the action of Eq. (3.23) can be easily transformed to Langevin equation for directed percolation 



d<i>(k,b;Y) ac 



dY 



2tt 



/ d 2 k'K(k,k')^(k',b;Y) 



2na 2 

-jj^*{k, b- Y)*(k, b; Y) + C(fc, b- Y) 



with 



and 



< C(k,b;Y) >= 



47T a 



< C(fc, b; Y)((k', b'; Y) >= $(fc, b; Y) 8 {2) (b - 9) 5^ (k - k 1 ) 8{Y - Y') 



N, 



(3.24) 



(3.25) 



(3.26) 



The Langevin equation is the one of many ways to describe a diffusion process and the considerable progress has 
been achieved in this approach (sec 



M, E3, 23, 37 ). 



However, we prefer a different way for description of the BFKL Pomeron interactions, which will also lead to 
diffusion equation: the so called generating functional approach. The advantage of the generating functional approach 
is its transparent relation to the partonic wave function of the fast hadron (dipole) . In this approach we see in the 
most explicit way our main theoretical problem: the BFKL Pomeron calculus provides the amplitude that satisfies 
the ^-channel unitarity while the s-channel unitarity is still a problem in the BFKL Pomeron calculus. However, the 
probabilistic interpretation in the framework of the generating functional leads to the correctly normalized partonic 
wave function which takes into account the main properties of the s-channel unitarity as well. 



3.3 Generating functional: general approach 



In this subsection we discuss the main equations of the BFKL Pomeron Calculus in the formalism of the generating 
functional, which we consider as the most appropriate technique for the probabilistic interpretation of this approach 
to high energy scattering in QCD. 

To begin with let us write down the definition of the generating functional ]17| 

Z(Y- Y ; [u]) = (3.27) 

~ n 

= ^ I P n (Y - Y Q ; xi,yi 1 ...;x i ,y i ;...- 1 x n ,y n ) Y[u(xi,yi)d 2 Xid 2 yi 

n=l i=l 

where u{xi,y{) = Ui is an arbitrary function of Xi and yi. The coordinates (xi,yi) describe the colourless pair of 
gluons or a dipole. P n is a probability density to find n dipoles with the size Xi — yi and impact parameter (a;, +yi)/2. 



It follows directly from the physical meaning of P n and the definition in Eq. (3.27) that the functional obeys the 
condition 



Z(Y - Y Q ; [u=l]) = 1. (3.28) 

The physical meaning of this equation is that the sum over all probabilities is equal to unity. 

Introducing vertices for the dipole process: 1 — ► 2 (Vi_»2(:r, y — > x\, yi+X2, 2/2)) and 2 — > 1 (V 2 ^i(xi, yi+x 2 , 1/2 - * 
x,y)) we can write a typical birth-death equation in the form 

& Pn\Y] . . . ; Xi, . . . 5 x n , 2/t2-J _ 



ay 



(3.29) 



E ^2 0(P„-i(F;.. )) (3-30) 

i 

E^ 2 ^i ( J P„+i(y;...;2; i ,y l ;...;x fe ,y fc ;...;a; n ,y n ) - P n (Y; . . . ; yv, . . . ;x k ,y k ; . ■ . ; x n ,y n )) (3.31) 



i>k 



Eq. ( 3.29 ) is the typical Markov chain and the fact that we have the correct normalized partonic wave function is 
written in Eq. Q3.29 ) by introducing for each microscopic (dipole) process two terms (see Eq. ( 3.30D and Eq. (3.31)): 
the emission of dipoles (positive birth term) and their recombination (negative death term). Multiplying this equation 
by the product Yii=i u i an< ^ integrating over all xi and yi, we obtain the following linear equation for the generating 
functional: 

dZ (Y-Y Q ;[u}) 



dY 



with 



= xM z (y-y ;[«D 



d qd qid q 2 [Vi 2 (q -> qi + 92) (-u(q) + u(q{) u(q 2 ) ) 



-V2-.1 (qi + qi -> q) (u(q 1 )u{q 2 ) - u(q)) 



Su(q) 



2 Su{qi)Su(q 2 ) 



(3.32) 

(3.33) 
(3.34) 



These equations describe the process of splitting of one Pomeron into two Pomerons and merging of two Pomerons 
into one Pomeron. We will prove below, that this set of equations is equivalent to the BFKL Pomeron calculus, 
given by the path integral of Eq. ( 2.1 3| ) . We use notations qi for (xi,yi) and d 4 qi for d 2 Xid 2 yi where Xi and yi are 
positions of quark (antiquark) of the colourless dipole. 



Eq. (3.32) is a typical diffusion equation or Fokker-Planck equation pM, with the diffusion coefficient which 



depends on u. This is the master equation of our approach, and the goal of this paper is to find the correspondence 
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between this equation and the BFKL Pomeron Calculus. In spite of the fact that this is a functional equation we 
intuitively feel that this equation could be useful since we can develop a direct method for its solution and, on 
the other hand, there exist many studies of such an equation in the literature ( see for example fl3|| ) as well as 
some physical realizations in statistical physics. The intimate relation between the Fokker-Planck equation, and the 
high energy asymptotic was first pointed out by Weigert Q in JIMWLK approach Q , and has been discussed in 

The scattering amplitude can be defined as a functional [p9l p"9| 

N(Y;M) = -£ J (~l) n ln(xi,y 1 ;...;x n ,y n ;Y )]J—Z(Y,[u t })\ Uz=1 d 2 x l d 2 y l (3.35) 



E( -1 )" / ln(xi,yi; ■ ■ .;x n ,y n ;Y ) p(xi, yr, ■ ■ ■ ; x n , y„; Y - Y ) tfxi cfyt 

n=l J 



The physical meaning of functions 7 n is the imaginary part of the amplitude of interaction of n-dipoles with 
the target at low energies. All these functions should be taken from the non-pertubative QCD input. However, in 
[EM 191 E0| it was shown that we can introduce the amplitude of interaction of n-dipoles 7„(xi, y\ \ . . . ; x n , y n ; Y) at 



high energies (large values of rapidity Y) and Eq. (3.27), Eq. (3.32) and Eq. (3.35) can be rewritten as a chain set of 
equation for j n (xi,yi; . . . ; x n , y„; Y). The equation has the form 5 



BY 



l] '"' ! ''' ' " qn) - = 2 E / A' A Vi 2 fe; q, q') 7 „ (...?'...)- E / <*Vi d ^2 9i, &) 7» (...,%.••) 

i— 1 z—l 

n— 1 „ n „ 

~ E / AcfVVi— 2(9*; q, g')7«+i (■■•<?• ■•</)- E / A^^i (ft, q ; q) 7n-i (%•••%••■ q) 



i=l 



i>3 



n „ n „ 

2 E / AA'^-i {q,qi;q') in-i (•••%• ■•<?) + E / A^-i 7« (•••&•••?*•••) 

(3.36) 



Comparing this equation for 71 = and 72 = T^ 2 ) one can see that 
Ki^2 = T^— r i^2 = —K(x,y;z) ; 

Z 7T Z 7T 

2a! a s / r 2 ^i(xi,yi + a; 2 ,2/ 2 -> 



7r 2 2 7r \ (a; — y) 4 

d 2 a; d 2 y 



(3.37) 
(3.38) 



(x - y) 4 



IW (a* , yi + x 2 , y 2 - y) (>> (a* - x)<5^ ( w - y) + c^ 2 > (x 2 - x)^ (y 2 - y)) 



with r 2 _^i is given by Eq. (2.35) 



3.4 A toy model: Pomeron interaction and probabilistic interpretation 

In this section we consider the simple toy model in which the probabilities to find n-dipoles being independent on 
dipoles sizes JIt], |l8|, |2(| In this model the master equation ( 3.32 ) has a simple form 



(3.39) 



5 This equation is Eq. (2.19) in [ pp| but, hopefully, without misprints, part of which has been noticed in ^?|. 



Eq. ( 3.39 ) generates: the Pomeron splitting Gp^2P = T(l — > 2); Pomerons merging G 2 p^p = |r(2 1) and 
also the two Pomerons scattering C?2P^2P = ^T(2 — ► 1). It is easy to see that by neglecting u 2 d 2 Z/du 2 term in 
Eq. ( 3.39 ) we cannot provide a correct sign for Pomerons merging G2P^p- It worth mentioning that we need this 
term to provide the conservation of the total probability on each level of rapidity. 



The description given by Eq. ( 3.39 ) is equivalent to the path integral of Eq. (3.8). To see this we need to notice 

(3.40) 



that the general solution of Eq. (3.39) has a form 

Z(Y;u) = e H{u ^ Y - Y ^Z(Y ;u) 
with operator H defined as 



and 



h{u) = _r(i-2)«(i-t0^+ ir(2-i) u (i-t0^ 



Z(Y ;u) = e^'"- 1 ' 



Introducing operators of creation (a + ) and annihilation (a) 



d_ 

du ' 



a + = u that satisfy [a, a + ] = 1 at fixed Y 



one can see that operator Ti. has the form 



n 



T(l -> 2) a+ (1 - a+) a + -r(2 -> 1) a+ (1 - a+) a 2 



and the initial state at F = Yq is defined as 



with the vacuum defined by a|0 >= 0. 



\Y > = e -^( a+ ~i)|0 > 



We need to discretize the development operator of Eq. (3.40) with TL given by Eq. (3.44), namely, 



JY 



e H(Y-Y ) = e HAY 



HAY 



= 11 (1 + HAY) 

3 = 1 



and introduce coherent states [B9| for a certain intermediate rapidity Yj = Yq + jAY as 



\<l>3 > 



|0> 



(3.41) 
(3.42) 

(3.43) 

(3.44) 

(3.45) 



(3.46) 



(3.47) 



where <pj are arbitrary complex numbers. The initial state of Eq. (3.45) can be written as 



\MYo) >= Vtr > ■ (3.48) 
The unit operator in terms of the coherent states can be expressed as 

1= (7 d</>jd * j e -^ r + ^ + *A |^ ><^| (3.49) 



We want to calculate matrix element of some operator A between states of initial Yq and final Y rapidity 
< F|j4|Yq >. This can be written as 



-i — > 



< Y\A\Y >=< Y\A 



(3.50) 



here we denote \Y >= \4>y >■ Next we use the development operator given in Eq. ( j3.46|) to find < Y\Yo >. We split 
the rapidity Y — Yq to N intervals and insert the development Eq. ( |3.46 ) and unit Eq. ( [3.49| ) operator between the 
states of intermediate rapidity 



JV 



We look at 



<Y\ JJ (1 + H AY) \Y > 

3=1 



< (1 + HAY) fa > = exp + cj>* +1 + - <j>* +l - <t> 3 + <f^ +1 4>j} 

x (i + w(^ +lJ 3 -)Ay) 

= exp - + ct> j+1 - (1 + H(<t>*+i,<l>3) AY ) 

= exp {-</>* +1 Oj + i - <£,■) + - <pj) exp Ay ) 



Now we redefine an arbitrary function <f)j as 

= -< 



and rewrite Eq. ( 3.52 ) in terms of <&j and 

e -^ + i(^+i-^)+^+i-^ exp {H(<j>* j+1 ,(j>j) AY} 
= exp {$ J+ i(4>+ +1 - $+) + H(l - -*+)} 



exp 



AY 



Summing over all rapidity intervals we have 

< Y\A\Y > ~ f[ { d<S>+d<S>jA{Y) 

3=0 J 

where A(Y) is the expectation value of the operator A at the final rapidity Y, and 

S = 



N 



*«^-*2 + „ (1 .» w ,_» ;) A y 



Ay 



In the continuous limit this becomes 



< y|A|y > = 



J X>$+X>$ A(y) e i 



with 



= [ dY ( $+-^-$ - r(i -> 2)$+$ + r(i -> 2)$+$ 2 + ir(2 -> i)$ +2 $ - ir(2 -> i)$+ 2 $ 



(3.51) 



(3.52) 



(3.53) 



(3.54) 



(3.55) 



(3.56) 



(3.57) 



(3.58) 



Eq. (3.58) is the action of Eq. (EDI) for T(l — ► 2) = \Y(2 — > 1). In the toy-model the difference between these 



two vertices is the normalization problem of functions <i> + and In our approach they are normalized in the way 



-i n 



which allows us to treat them as probabilities (see Eq. ( |2.17 ) and Eq. ( 2.18 )). However, Eq. ( 3.5S ) includes the 
new interaction: a transition of two Pomerons to two Pomerons. The sign is such that this interaction provides the 
stability of the potential energy. Indeed this term is responsible for the increase of the potential energy at large 
values of both $ + and <E>. 

Comparing Eq. ( 3.58| ) with Eq. (3.1) one can see that we built the partition function and the thermodynamic 
potential using the generating functional. It means that our Eq. (3.32) is equivalent to statistical description of the 
system of dipoles. 

Eq. (3.39) is the diffusion with the u dependence in diffusion coefficient. In terms of the Langevin equation 
Eq. ( |3.58 ) generates a noise term of the < >oc $(1 — $) type. 

To out taste Eq. ( 3.3S ) is simpler than the Langevin equation of Eq. ( |3.13|) and it will be easily generalized for 
the case of QCD. For u < 1 the diffusion coefficient is positive and the equation has a reasonable solution. If u > 1, 
the sign of this coefficient changes and the equation gives a solution which increases with Y and Z(Y) cannot be 
treated as the generating function for the probabilities to find n dipoles (Pomerons) (s ee |To| , |l2| , pl| for details) . The 
same features one can notice in the asymptotic solution that is the solution to Eq. ( 3.39 ) with the l.h.c. equal to 
zero. It is easy to see that this solution has the form 



Z(Y — > oo; u) — 



1 - 



with k = 2 



r(i - 2) 
r(2 - 1) 



2 7V 2 



> 1 



(3.59) 



For negative k this solution leads to Z > 1 for u < 1. This shows that we cannot give a probabilistic interpretation 
for such a solution. 



4. A practical way to find solution: Monte Carlo simulation 



In this section we consider the BFKL Pomeron calculus in the form of Eq. (3.22) which leads to the simplest approach 
in the framework of the generating functional technique. This approach is not only the simplest one but also it is 
free from all troubles related to the negative contribution for the process of 2P — > P transition. Repeating procedure 
discussed in section 3.2 one can see that Eq. ( 3.22 ) for action leads to the following equation for the generating 
functional 

dZ (Y-Y ;[u(b,k)]) 



dY 



with 



xm 



d 2 bd 2 k 



as 
2tt 



X [u(b,k)} Z (Y-Y 0] [u(b,k)\) 



d 2 k' K{k, k')u(b, k') - , + u(b, k) u(b, k) 



4:Tr 2 as\ as 



5u(b, k' 
1 S 2 



Su(b, k) 



27T (u(b,k)u(b,k)-u(b,k)) 2Hbjk)Su(b k) 



For Z(Y, [u{b 7 k)}) defined as 



Z(Y - Y 0] [u(h,ki)}) = 
^ / Pn(Y-Y ;b 1 ,k 1 ;...;bi,k i ;...;b n ,k n )Y\_u{bi,k i )d 2 



(4.1) 

(4.2) 
(4.3) 

(4.4) 



where u{bi, k{) are arbitrary functions 
Eq 



and Eq. (4.3) show that the 2P —* IP transition can be written as the two dipole to one dipole merging 
with a positive probability. Therefore, Eq. (4.1) has a very simple probabilistic interpretation which can be written 
as the following Markov chain: 

d Pn\Y] • • • 5 b$ = b 7 ki = fej . . . j b n — 6, k n = &) 



dY 



(4.5) 



-i f\ 



— [Pn-i(Y] ...\bi = b,ki = k;...; & n _i, fc n -i) - 

i 

- J d 2 k' K(k, k') P n (Y; . . . ; h = b, k t = k'; . . . ; x n = b, k n = k)^j 
4 7r as \ Q-s 



(4.6) 



) ff (- P «+ 1 C^> • • • ' h = b, kj = k; . . . ; h = b, h = k; . . . ; b n = b, k n = k)- 



i>j 

- P n (Y; . . . ;bj = b,kj = k; . . .h = b, ki = k; . . . ;b n = b, k n =k;)); (4.7) 

This set of equations can be solved numerically and it gives a practical way to discuss the influence of the 
Pomcron loops on the solution for the scattering amplitude at high energies. 



5. Conclusions 



We demonstrate in this paper that the BFKL Pomeron Calculus in the kinematic region given by Eq. (1.3) has 
two equivalent descriptions: (i) one is the generating functional which gives a clear probabilistic interpretation of 
the processes of high energy scattering and provides also a Hamiltonian-like description of the system of interacting 



dipoles; (ii) the second is the Langevin equation with a noise term of Eq. (3/7) which is rather complicated. We 



show that at high energies the Langevin equation with noise of Eq. (3.7) can be reduced to the Langevin equation 
for directed percolation in the momentum representation if the impact parameter is large, namely, 6 > 1/i where k 
is the transverse momentum of a dipolc. Unfortunately, this simplified form of Langevin equation is not applicable 
for summation of Pomeron loops, where one integrates over all possible values of impact parameter. 

In other words, the BFKL Pomeron Calculus can be considered as an alternative description of the statistical 
system of dipoles with different kinds of interactions between them. The same conclusion was drawn in |5q ] using a 
similar technique. We show that the BFKL Pomeron calculus with two vertices: splitting of one Pomeron into two 
Pomerons(P — * 2P) and merging of two Pomerons into one Pomeron (2P — > P), can be described as a system of 
colourless dipoles with two processes: the decay of one dipole into two dipoles and the merging of two dipoles into 
one dipole. 

It is shown that the question about negative amplitude does not arise if we treat the system of dipoles in 



the momentum representation. Markov chain for this system is written in the paper (see Eq. (4.5)) and can be 
considered as a practical way to find a solution in accessible range of energies. However, as we have stressed a 
number of assumptions we have made on the way to simplify the noise term as well as to argue the probabilistic 
interpretation of the Pomeron calculus. Therefore, both these approaches can be considered as a QCD motivated 
models. In the case of the probabilistic interpretation we can even expect that this approach will not work in the 
form, that has been suggested in this paper, for the full (2P — ► P) vertex (see p7fl). 

Being elegant and beautiful the BFKL Pomeron Calculus has a clear disadvantage: it lacks theoretical ideas 
what kind of Pomcron interactions wc should take into account and why. Of course, Feynman diagrams in leading 
ln(l/cc) approximation of perturbative QCD allow us, in principle, to calculate all possible Pomeron interactions but, 
practically, it is very hard job. Even if wc calculate these vertices we need to understand what set of vertices we should 
take into account for the calculation of the scattering amplitude. This is the reason why we need to develop a more 
general formalism. Fortunately, such a formalism has been built and it is known under the abbreviation JIMWLK- 
Balitsky approac h ]4l[ In this approach we are able to calculate all vertices for Pomeron interactions as it was 

demonstrated in fig] and it solves the first part of the problem: determination of all possible Pomeron interactions. 
However, we need to understand what vertices we should take into account for calculation of the scattering amplitude. 
We hope that a further progress in going beyond of the BFKL Pomeron Calculus (see [fh| 0) will lead to such a 
development of the BFKL Pomeron Calculus with a consistent theoretical approach. Hopefully this approach will 
be simpler than Lipatov effective action E§] which is not easier to solve than the full QCD Lagrangian. 



It is well known that the mean field approach to our problem, which includes only one dipole to two dipoles 
decay in the master equation (see Eq. ( [3.32 )) has been studied quite well both analytically plj and numerically f33|| . 
We firmly believe that the probabilistic interpretation in spite of being a QCD motivated model, gives a practical 
method for creating a Monte Carlo code in spirit of the approach suggested in Q. This code will allow us to find 
a numerical solution to the problem and to consider inclusive observables. This extension is very desirable since the 
most experimental data exist for these observables. 



Acknowledgments 

We want to thank Asher Gotsman, Larry McLerran, Dima Kharzeev, Alex Kovner, Misha Lublinsky and Uri Maor 
for very useful discussions on the subject of this paper. Our special thanks go to Jose Guilherme Milhano, who draw 
our attention to |3^, |35| and the discussions with whom on the subject of this paper were very useful and instructive. 

This research was supported in part by the Israel Science Foundation, founded by the Israeli Academy of Science 
and Humanities and by BSF grant # 20004019. 

A. Calculation of G (x 1 , x 2 \x' 1 , x' 2 ) 

The solution of the BFKL equation is given by [^6| 

, / f +oc {v 2 +n 2 /A)dv G ufi {x 1 ,x 2 ,x' 1 ,x' 2 ) 
G{x 1 ,x 2 ;x 1 ,x 2 \uj) = > / TlWm — 7 ^ ( AJ ) 



where G Vfl (xi, X2,Xi,x 2 ) is the Mellin transform of Eq. ( |2.q ) and Xi are two-dimensional vectors in complex 
coordinates 



,y — ^l^x t^i^y 



(A.2) 



The function u>(i>,n) is the eigen value of the BFKL equation given by Eq. (|2.8|). 

The four-point Green function is presented in terms of the hypergeometric functions p6| , ^| 

G UfJ ,(xi,X2,x 1 ,x 2 ) = C\x h x* h F(h, h, 2h; x)F(h, h, 2h; x*) 
+C 2 x 1 - h x* 1 ~~ h F(l -h,l-h,2-2h; x)F(l -h,l-h,2-2h;x*) (A.3) 

with h = h + iv + t- , h — | + iv — ^, and x is the anharmonic ratio 

X12XV2' , 

x = (A. 4) 

X 1V X 2 2' 



Coefficients C\ and C2 are given by 26 



2tt 2 ~' 2ir 2 



with 



,o4^ r (-^+(l + W)/2)r(^+|n|/2) 



bn '"~ w 2 '"% + (i + H)/2)r(- w +H/2) (A ' 6) 
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As we have discussed, the high energy asymptotic behaviour stems from n = term in Eq. (A.l). The initial 
condition for Eq. ( |2.6|) at Y = Yq is given by the following expression 



Go(xi,X2;x' 1 ,x' 2 ) = (A. 7) 

dui f +OQ v 2 dv G v ^{xi,X2,x' 1 ,x' 2 ) 1 



27riJ_ 00 (^ 2 + l/4) 2 lu -w{v,n = Q) ^ J-oo {v 2 + l/A) 2 
This integral can be taken by closing contour of integration over singularities of the integrand 
The function 



x 2 ,x' 1 ,x' 2 ) 



1/4)2 



(A., 



has two poles at | and —&. The four-point Green function of Eq. (A. 2) consists of two terms, one with Ci(xx*)i +W 
and the other one C 2 {xx*) 3 ~ lv . For small |x| this terms should be integrated closing contour in upper and lower 
semi-planes respectively. The resulting contour in the lower semi-plane runs anticlockwise and thus the value of the 
contour integral enters with a minus sign. 

The terms could be expanded in the vicinity of their poles. Let us consider the first term. We expand the 
function G\ in the vicinity of | 

The hypergeometric function can be written as a sum 

T, f r s. r(c) ^ r(o + n)T(b + n) x n in , 

F(a, b,c;x) = l+ y v ; v ^— A.10 

r W r W^ r(c + n) n! 

In the case of -F(/i, h, 2/i; x) the singularity of F function at | can factorized out the sum 

F(fe,/ t ,2fe;x) = l+ N 



r(fc)r(/o^ r(2/i + n) ni 



1 ^ x™ 1 ira(l - x) 



Vi_ = l-i '^-^ (A.n) 



2r(/i)^n 2r(u/+§) 
At this stage the first term of the integrand of Eq. ( |A.7| ) can be written as 



(A.12) 



w (w + |) 3 (w- 1) 2 2 r(-w + |)r(i + ; V 2r(i + w)yV 2r (| + ^) 



It is clearly seen that the term of zero order in ln\l — x\ has a third order pole at v — ^; the term of first order 
in ln\l — x\ has a second order pole at v = ~, and, the term of second order in ln\l — x\ has a simple pole at v — ^. 
The contributions of those terms in the contour integral are found to be 



.2 



-2 + /n(ra*)(-4 + ;n(xa;*))] (A.13) 



J- [-2 + 7 + In(a;3*)]Jn[(l-a;)(l-a:*)] (A.14) 



-iuJ 



ln(l - x)ln{l - x*) 



respectively. 

In a similar way we may expand the second of the integrand in vicinity of its pole at v = — ^ namely 



(A.15) 



-~2 



(xx 



,.- , 1 _lln(l-x) 



2r(|-ii/) 



lZra(l-a;*)' 
~2r(±-u/), 



(A.16) 



The integration is performed on the lower semicircle and results in overall minus sign of the integral. The contributions 



corresponding to Eq. ( A. 13 ), Eq. ( A. 14 ) and Eq. ( A.15 ) are 



[-2 + ln{xx*){-A + ln(xx*))] 



2uj 



-2 + 7 + ln(xx*)]ln[{l - x)(l - x*)] 



+ —ln(l - x)ln(l -~ x*) 



respectively. 

Comparing the contributions we note that that of zero and second order 
in ln[(l — x)(l — x*)] are exactly canceled out, and we are left with 



For small led this can written as 



2 — [-2 + 7 + ln(xx*)]ln[(l - x)(l - x*)] 
2lu 



4tt 2 

In \x\ lull — x\ 

u> 



4tt 2 , 1 , , 
In- — r In 1 — x\ 

LU \X\ 



Going back to complex vector representation of x and rewriting Eq. (A. 23) as 

47T 2 

Go(xi,X2\x 1 ,x 2 ) = H In 



4tt 2 


In 




In 


XWXll' 












X\2>X V 1 




X12XV2- 



(A.17) 
(A.18) 
(A.19) 



(A.20) 

(A.21) 
(A.22) 

(A.23) 



we see that we reproduce the result of |26|]. Therefore, we demonstrated that Eq. ( A.23 ) gives a correct initial 
condition for searching the scattering amplitude at high energies restricting ourselves by the one term in Eq. ( |2.6| ) 
with n = 0. 



B. The path integral formalism for the generating functional 

We want to develop a path integral formalism similar to that we found for the toy model, but where the probabilities 
to find n dipoles depend of dipole sizes. In our notation we denote by Latin index rapidity interval, and Greek indices 
relate to a size of dipole. 

As in Chapter 3.3 we introduce the creation and annihilation operators 



&{q)= S^(q} ^(q)=u(q) (B.l) 

with commutation relations [a(q),a'(q')] = S(q — q') at fixed Y. The expression for the coherent states in this case 
takes form of 



\4>{q) >=e^ a |0> (B.2) 

with 

a{q')\ct>{q) >= <t>(q)\<t>(q) > 6(q - q') (B.3) 

First we consider discrete dipole sizes q a = (L/N)a, where L is a maximal possible dipole size, TV a number of 
intervals of L discretization, and a is an integer number running from to N. In this case the commutation relations 
become [a(q a ), a^{qp)] = 5 a ,p- 

The unit operator can be written in terms of the coherent states 

I = TT / MgjMgg) e -^(fa)»(fa)+^(fe)+^(fa)|0( go ) >< ^ qa )\ (B.4) 

a J 

The operator TL defined in Eq. ( 3.33| ) can be written as 



H[a\a}= - ^2^2^2[Vi^2(qp -> q~, + qx){-^{q p ) + (q^tf (qx)}a(qp) 
7 a 

- V 2 ^i(qj + q\^ qp){a J {q-,)a) (q x ) - a T (9 /3 )}ia(g 7 )a(g A )] (B.5) 
Following the logic of Chapter 3.3 we consider a matrix element 

jn < &+i(?«OI j (1 + HAY) 1 1 4>j(Qa) > j (B.6) 
First, look at the second term of the Hamiltonian 

\ II < <t>i+i(9<*')\ \J2J2J2^- Vl ^ 2 ^ *r + 9A)a t (g 7 )a t (gA)a(^)] ^ ]J |0,(g Q ) > I 

la' J/37A la J 

la' J/37A la J 
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In Eq. (B.7) we used the property of the coherent states given by Eq. ( |B.3| ). 

In the continuous limit S a ^ is replaced by 5(q a — qp), and Y[ a d(f>* (q a )d(f>(q a ) by functional integration J T>(p*'D<f). 

From here we see that rest of the calculations is similar to that of Chapter 3.3 and we end up with the expression 
for a matrix element of an operator A between states of initial Yq and final rapidity Y 



where 



S 



< Y\A\Y >~ J V<S> + V$A{Y)e s 
J (/ 3 + (9)^*(«)*f + W( 1 -*>-* + )W 



with the Hamiltonian given by 



(B. 



(B.9) 



H = 



d 4 q d 4 q 1 d 4 q 2 [V 1 ^ 2 (q -> q x + q 2 ){-<&{qo) + $(<Zi) + - $(<7i)$(<?2)}$ + 0?o) 



+ V 2 ^i{qi +q 2 -> «b){-*(«l)*(«2) + + - $(<?o)}2* + (9i)* + (<72)] 



(B.10) 
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